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Abstract 

We combine the results of CDF and D0 searches for the rare decays B® — > [i^ [i~ and 
Bj — > fi + fi~ . The experiments use 364 pb~ x and 300 pb _1 of data respectively. The 
limits on the branching ratios are obtained by normalizing the estimated sensitivity 
to the decay B + — > J/ipK + taking into account the fragmentation ratios f u /f s (d)- 
The combined results exclude branching ratios of BR(5j — > > 1.5 x 10~ 7 and 

BR(B® — > fi + fi~) > 4.0 x 10~ 8 at 95% confidence level. These are the most stringent 
limits on these decays at the present time. 
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1 Introduction 



The CDF and D0 experiments have previously reported on searches for the rare decay B® — > 
A i+ A t ~ |H E|- CDF also directly searched for the decay B® — > These decays are highly 

suppressed in the Standard Model of particle physics with branching ratios of BR(5° — > 
= 3.5 x and BR(B® — > ,u + /i~) = 1.0 x 10~ 10 However, decays of this type can 
be significantly enhanced in many scenarios beyond the Standard Model [4]. A combination 
of results leads to more stringent limits and is of considerable interest in exploring the phase 
space of the models where strong enhancements for B® — > /i + jjT or B® — ► are predicted. 

In this note we report on a combination of limits in the B® — > and B^ — > n + fi~ decay 

channels. 

CDF and D0 use similar methodologies to search for the B° s — * decay. CDF applied 

the same methods to directly search for the decay B® — > Each experiment looks for two 

oppositely charged muons in the B° s and B% mass range using dedicated triggers. CDF divides 
their dataset into two channels: the " Central" channel consists of muon pairs reconstructed in 
the pseudorapidity region, \r)\ < 0.6, and the " Central- Extended" channel consists of dimuon 
events where one muon is reconstructed in the central region and the second muon in the 
extended muon system, 0.6 < \r]\ < 1.0. The two channels have different sensitivities, therefore 
the optimization is performed separately for each channel. The branching ratio is computed 
by normalizing the number of signal events to the number of reconstructed B + — ► J/ vf)K + [3] 
events. The branching ratio or limit is then calculated from the equation: 

ATobs total £ 

BR(5° - ^V) = -^%BT ■ ■ T ■ BR ( 5+ - J /^ + ) • BR W "> A*V), (1) 

where: 

• Ng^ is the number of observed B° s — > n + \i~ candidates; 

• «b s is the geometric and kinematic acceptance of the di-muon trigger for B° s — > fJ, + /J,~ 
decays; 

• e^ s tal is the total efficiency (including trigger, reconstruction and analysis requirements) 
for B® — > ji + \T events in the acceptance; 

• N B +, ccb+, and e^+ al are similarly defined for B + — > J/tpK + decays; 

• fu/ fs accounts for the different 6-quark fragmentation probabilities and is: 

(0.397 ± 0.010)/(0.107 ± 0.011) = 3.71 ± 0.41, where the anti-correlation between 
the uncertainties has been accounted for 0; 

• BR(5+ -> J/i)K + ) ■ BR(J/ip -> p+fT) = (1.00 ± 0.04) x lO" 3 ■ (5.88 ± 0.10) x 10~ 2 = 
(5.88 ± 0.26) x 10~ 5 are used . 

The experiments normalize to the decay mode B + — > J/ ipK + rather than to the B s decay 
B° s — > J/ip<f). Normalizing to the decay B + — > J/ipK + is preferable since the mode has 
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higher statistics and the branching ratio and lifetime are well known from the measurements 
at CLEO and the asymmetric B factories. In addition, understanding the efficiency to detect 
— > J/ ip<p events is complicated by the presence of CP even and odd decay components which 
have different lifetimes. Finally normalizing to the mode B° s — > J/ip(j) does not eliminate the 
systematic uncertainty from the ratio f u / f s since current measurements of the branching ratio 
of B° s — > J/ip(p are calculated using the fragmentation ratio. The expression for B^ — > 
is obtained by replacing B® with B® and the fragmentation ratio with f u j which is taken as 
unity. The limits on the branching ratio BR(5° — > are also calculated using the ratio 

fu/ fs [Z] determined from Tevatron data based on the CDF Run 1 analysis jH]. 

To increase the sensitivity both experiments perform an optimization over the primary 
discriminating variables, which are based on lifetime, compatibility between the momentum 
vector of the candidate B meson and the vector between the production and decay vertices, 
and isolation in a cone around the candidate meson, where isolation is defined as the ratio 
of sum of momenta of the candidate tracks divided by the sum of momenta the candidate 
tracks and other tracks in the cone. CDF optimizes for the best limit, and D0 for the best 
sensitivity, where sensitivity is defined as e Bs /(l + ^/Nbg)-, and Nbg is the number of expected 
background events. 

We use a Bayesian integration method to calculate the combined limits jO]. The method 
takes into account correlated and uncorrelated systematic uncertainties between the two ex- 
periments, and between the two search channels of the CDF analysis. 

The key elements in calculating the limit are estimating central values and uncertainties 
on the efficiency and acceptance for the trigger and reconstruction, and estimating the back- 
ground. These elements are described in the next section, followed by a section on the limit 
calculation and a final section where the results are summarized. 

2 Acceptance, Efficiency, Backgrounds and Associated 
Uncertainties 

Both experiments evaluate the acceptance and efficiency of their respective triggers, recon- 
struction code and discriminating variables. By normalizing the limit to the measured decay 
B + — ► J/%l)K + most of the systematic uncertainties of these estimates cancel in the ratio. 
The sources of systematic uncertainties in the CDF and D0 analyses are discussed in detail 
in references Systematic uncertainties on the estimates of the acceptance and efficiency 

are uncorrelated since they are based on estimates of trigger and reconstruction efficiencies for 
different detectors. Systematic uncertainties associated with the production of B mesons are 
possibly correlated since similar event generation programs are used in both analysis. These 
systematic uncertainties affect the estimates of the efficiencies of the primary discriminating 
variables. However, the variables that are used are qualitatively different. For instance, the 
strongest discriminator in each analysis is the lifetime variable. The CDF analysis uses the 
proper decay time in 3D while D0 selects events based on the 2D transverse decay length 
significance distribution. The discriminating variable used by D0 is correlated with mo- 
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mentum where the CDF variable is uncorrelated. Therefore the uncertainties on f u /f s and 
BR(5 + — ► J/i/jK + ) ■ BR(J/^ — > /i + yU~) are 100% correlated between the two experiments 
and all other uncertainties are treated as uncorrelated. 

The experiments evaluate the background using sideband events. The uncertainties on the 
backgrounds are dominated by the statistical uncertainties. The acceptance, efficiency, and 
background numbers are summarized in Tabled 





CDF: Central 


CDF: Central-Extended 


D0 


Luminosity 


364ptT i 


336pb 1 


300pb~ x 


total 

,a B+ -e B+ . 

V „ total ) 


0.852 ±0.084 (±9.9%) 


0.485 ±0.048 (±9.9%) 


0.247 ±0.019 (±7.7%) 


ATOfcS 

iV B+ 


1785 ± 60 (± 3.4%) 


696 ± 39 (± 5.6%) 


906 ±41 (±5.0%) 


Uncor. Uncer. 


(±10.4%) 


(±11.3%) 


(± 9.2%) 


fu/fs 


3.71 ±0.41 (±11.0%) 


BR(B+ -► J/i)K+ 


(5.88 ± 0.26) x 10- 5 (±4.0%) 


Cor. Uncer. 


(±11.9%) 


Kjexpected 
ly back 


0.81 ±0.12 (±14.8%) 


0.66 ±0.13 (±19.7%) 


4.3 ±1.2 (±27.9%) 


ivrobs 
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ses (xlO 7 ) 


1.04 ±0.16 (±15.8%) 


1.52 ±0.25 (±16.4%) 


0.59 ±0.09 (±15.0%) 


ses (xlO 7 ) 


0.617 (CDF combined) 




Exp Limit 90% CL 


3.5 x 10- y 


5.6 x 10~ 7 


3.5 x 10~ Y 


Exp Limit 90% CL 


2.0 x 10" 7 (CDF combined) 





Table 1: A summary of the inputs used in equation^ to estimate BR(5° — > fi + fi~). The 
relative uncertainties are given parenthetically. The single-event-sensitivity, ses, to a given 
branching ratio, corresponding to N?^ 1 = 1, and the expected limit at 90% Confidence Level 

^ s 

(CL), under a hypothetical repetition of the experiments, are calculated using the inputs. 
The combined ses and expected limit for the CDF "Central" and "Central-Extended" search 
channels are also given. 



3 Limits 

Using the Bayesian integration method we calculate the combined limits. In the case of the 
D0 search the dimuon mass signal region covers both the and B%. The limit on the 
branching ratio that is extracted in one mode assumes that the branching ratio in the other 
mode is zero, which results in a conservative limit. This is the case in the framework of Minimal 
Flavor Violating (MFV) SUSY models, where the CKM matrix is the only source of flavor 
violation. In MFV SUSY the branching ratio for B^ — > /i + pT is expected to be suppressed 
relative to B® — > ji + by a factor of (V^/V^l 2 , making the contribution from B\ — > 
negligible. The 95% confidence level (CL) limits on the branching ratio BR(5° — > are 
reported in Tabled The combined limits are also calculated at 90% CL. 
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CDF: Central 


CDF: Central-Extended 


D0 


T ill m inrmitv 

1 ' LLllllllWOl Ij y 


364 pb" 1 


336 pb" 1 


300 pb" 1 


Limits 
at 95%CL 


BR(B° S -> < 2.0 x 10~ 7 


BR(B° S -> < 3.9 x 10~ 7 


Combined 
Limits 


BR(B° S -> At+//-) < 1.5 x HT 7 at 95% CL 
BR(5° -> < 1.2 x HT 7 at 90% CL 



Table 2: A summary of the limits on the branching ratio BR(5° — > ). 



For the branching ratio BR(i?° — > the factor f u /fd is taken to be unity. The 

efficiency for the S° channel is estimated to be 8% lower than in the 5° channel for the 
D0 search, and the uncorrelated uncertainty is 10.2%. The CDF efficiency and errors are the 
same as for the B® —>■ fi + n~ channel. The limits on the branching ratio BR(B d — ► are 
reported in Table El 





CDF: Central 


CDF: Central-Extended 


D0 


Luminosity 


364 pb~ x 


336 pb" 1 


300 pb" 1 


Limits 
at 95%CL 


BR(B° d -»■ < 5.1 x 10~ 8 


BR(B° d -> ^/x-) < 11.1 x 10~ 8 


Combined 
Limits 


BR(B° d -> < 4.0 x 10~ 8 at 95% CL 
BR(52 -> < 3.2 x 10" 8 at 90% CL 



Table 3: A summary of the limits on the branching ratio BR(B d — > ). 



The limits on the branching ratios BR(B® — > and BR(B d — > can be compared 

to the values expected in the SM. Currently the 95% CL limit on BR(5° — ► /i + /U~) is approx- 
imately a factor of 50 larger than the SM branching ratio while the limit on BR(B d — > a + u~) 
is approximately 400 times larger. 

The uncertainty on the ratio f u / f s is the largest systematic uncertainty in limit calculation. 
This error dominates the external sources of systematic uncertainties. In order to facilitate 
recalculation of the limit as the measured value of f u / f s is improved we give the limits factoring 
out this contribution in Table EJ 

The limit on branching ratio BR(B® — ► can also be calculated using an evaluation 

of the fragmentation fractions based only on Tevatron data. This average, f u /fs = (0.398 ± 
0.010)/(0.120 ± 0.021) = 3.32 ± 0.59, differs slightly from the world average value and is 
dominated by the comparison of the mixing probabilities between B® and B d mesons, Xs and 
Xd- Limits calculated using this value of f u / f s are reported in Table El 

The limits were checked using independent programs for Bayesian limit integration devel- 
oped separately by the CDF and D0 collaborations. In addition, the cutoff of the assumed 
prior distribution was varied and the distributions used for the uncertainties in the calculation 
were taken as Gaussian distributions with a cutoff and Gamma functions. The limits were 
found to vary by less than 5% under these tests. 
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CDF: Central 


CDF: Central-Extended 


D0 


T ,n in i n OQi"t""\r 
u LiliiiiiUol u y 


364 pb" 1 


336 pb~ x 


300 pb" 1 


Limits 
at 95%CL 


BR(5 S ° - fJL+fJL-yifJfs) < 5.1 x 10- 8 


BR(B° -> iM + iM-)/(f u /f.) 
< 10.1 x 10- 8 


Combined 
Limits 


BR( J B S ° -> fi + fi~)/(f u /fs) < 3.9 x 1(T 8 at 95% CL 
BR(B° S -> /iVi/iy/j < 3.1 x 10- 8 at 90% CL 


Table 4: A summary of the limits on the branching ratio BR(B® — > )/(f u /f s )- 




CDF: Central 


CDF: Central-Extended 


D0 


Luminosity 


364 pb" 1 


336 pb" 1 


300 pb" 1 


Limits 
at 95%CL 


BR( J B S ° -> < 1-9 x 1(T 7 


BR(B° S -> ^ V) < 3 - 7 x 10_7 


Combined 
Limits 


BR(B° -> < 1.4 x 10~ 7 at 95% CL 
BR(£° -> At+At") < 1-1 x 10~ 7 at 90% CL 



Table 5: A summary of the limits on the branching ratio BR(B® — > ) using the Tevatron 
determination of f u / f s . 



4 Conclusion 

We report on a combination of CDF and D0 limits on the branching ratios BR(£?° — > / u + 
and BR(5° — > fi + fi~). The limits are obtained using a relative normalization to the decay 
B + — > J/ipK + . The individual limits are combined using a Bayesian integration technique 
that takes into account correlated and uncorrelated systematic uncertainties between the two 
experiments. The limits are found to be robust under several tests. The combined limits are: 

7 at 95% CL 

7 at 90% CL 

8 at 95% CL 
8 at 90% CL 



BR(B° - 




> 


BR(B° S - 




> 


BR(B° d - 




> 


BR(B° d - 





) < 1.5 x 10- 
) < 1.2 x 10" 
) < 4.0 x 10" 
) < 3.2 x 10- 
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